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ABSTRACT 


The  attenuation  of  seismic  waves  in  the 
crust  and  upper  mantle  has  been  studied 
by  various  techniques  in  order  to  determine 
the  specific  quality  factor  Q.  Anelastic 
deficiencies  of  the  propagation  medium 
can  only  be  separated  after  corrections 
are  made  for  the  effects  of  instrumental 
distortion,  of  the  "receiver  crust"  and 
the  "source  crust",  the  lossless  propa¬ 
gation  through  the  mantle  and  the  source 
mechanism. 

Seismic  refraction  data  have  been  compiled 
for  a  number  of  "source"  and  "receiver" 
regions  in  Europe.  A  refined  velocity-depth 
distribution  with  two  velocity  reversals 
and  an  average  value  of  Qp^670  have  been 
obtained  for  the  crust  of  the  Bohemian 
Massif  in  Central  Europe.  The  Qp  of  the 
uppermost  mantle  was  found  to  bo  about 
250,  Crustal  transfer  functions  were 
calculated  to  test  models  of  crustal 
structure  derived  by  other  methods.  The 
radiation  patterns  for  two  types  of 
explosive  sources  have  been  determined. 
Their  shape  is  strongly  affected  by  the 
source  mechanism  and  crustal  structure. 

A  now  velocity-depth  structure  of  the 
upper  mantle  in  Europe  has  been  derived 
from  travel-time  observations  of  large 


explosions  and  earthquakes .  At  least  three 
velocity  reversals  have  been  found  down 
to  depths  of  about  220  km.  They  correlate 
well  with  the  distribution  of  earthquake 
focal  depths.  Surface-wave  dispersion 
measurements  are  also  in  agreement  with 
these  now  findings.  Amplitude  decay  data 
for  surface  waves  indicate  that  zones  of 
low  Qg  must  exist  at  depths  where  the 
low-velocity  channels  were  found. 

Finally  an  average  value  of  Q^«220  for 
the  deeper  mantle  has  been  obtained  based 
on  observations  of  two  intermediate  focus 
earthquakes.  A  comparison  with  corresponding 
(V  data  shows  that  Qn«*Qe,  i.e.  that  losses 
in  pure  compression  cannot  be  neglected. 
Alternatively  it  is  suggested  that  Q  may 
be  frequency-dependent,  at  least  over  an 
extended  frequency  band. 
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INTRODUCTION 


1*1  Previous  Results 

Attenuation  of  seismic  waves  has  been  measured 
for  many  years*  But  it  is  only  recently  that 
sufficient  precision  in  observation  has  been 
achieved  and  that  analytical  techniques  have 
been  developed  which  make  possible  the  inversion 
of  the  experimental  body  and  surface  wave  data 
(ref.  KNOPOFF,  1964}  ANDERSON  an.  KOVACH,  1964} 
ANDERSON  and  ARCHAMBEAU,  1964}  BEN-MENAHEM, 

SMITH  and  TENG,  1965)*  The  resultant  distri¬ 
bution  of  attenuation  versus  depth  for  com- 
pressional  and  shear  waves,  usually  described 
by  the  dimensionless  quality  factors  Qp  and 
Qs,  is  an  important  and  valuable  source  of 
additional  information  with  regard  to  the 
pressure,  temperature,  state  and  composition 
in  the  earth's  Interior.  It  supplements  the 
velocity  data  and  provides  a  measure  of 
anelasticity,  i.e.  it  permits  a  description 
of  the  deviation  from  perfectly  elastic  wave 
propagation. 

Results  of  seismological  investigations  over 
the  past  decade  have  demonstrated  beyond  doubt 
that  the  structure  of  the  earth's  crust  and 
upper  mantle  is  much  more  complex  than  has 
‘been  assumed  until  recently.  It  is  no  longer 
permissible  to  approximate  the  uppermost  300 
to  400  km  of  the  earth  by  a  model  consisting 
of  a  few  isotropic,  homogeneous,  lossless  flat 
layers,  as  is  for  instance  still  done  in  the 
interpretation  of  a  number  of  geophysical  data. 


Phase  velocity  measurements  of  seismic  surface 
waves  have  shown  that  beneath  a  highly 
differentiated  crust  a  zone  of  low  shear- 
velocity  (Vg)  must  exist  in  the  upper  mantle 
at  depths  between  100  and  200  km  (see  Fig*  1 , 
right) •  This  zone  which  is  associated  with 
temperatures  approaching  the  melting  point 
of  the  material,  probably  has  a  higher  elec¬ 
trical  conductivity  than  the  regions  above  and 
below  it*  Its  upper  boundary  which  in  the  Alps 
and  in  their  northern  foreland  is  found  at  a 
depth  of  about  80  km  (ref.  KNOPOFF,  MUELLER 
and  PILANT,  1966}  SCHNEIDER,  MUELLER  and 
KNOPOFF,  1966)  must  be  identical  with  the 
discontinuity  which  supposedly  separates  the 
resistosphere  from  the  conductosphere*  The 
seismic  velocities  in  this  asthenosphere 
channel  which  extends  to  a  depth  of  about 
200  to  220  km  are  lowered  by  about  ten  per 
cent*  It  will  be  shown  subsequently  that  this 
simple  picture  is  much  too  crude  and  requires 
considerable  modifications* 

In  addition  to  these  Observations  evidence 
has  been  presented  that  a  distinct  reversal 
of  compressional  velocity  (Vp)  must  occur  at 
depths  between  about  7  and  12  km  within  the 
continental  crust  (see  Fig*  1,  left)*  The 
Interpretation  of  strong  second  arrivals  on 
seismic  refraction  profiles  as  well  as  of 
large  normal-incidence  reflections  requires 
the  presence  of  a  well-developed  low-velocity 
zone  in  the  sialic  part  of  the  crust  overlying 
a  layer  whose  velocity  is  markedly  greater  than 
the  velocity  in  the  channel  and  at  least  0.2  km/ 
seo  higher  than  that  above  the  channel  (ref* 
MUELLER  and  LANDISMAN,  1966)  . 


Fig,  It  Sumary  of  velocity-depth  distributions 
in  the  crust  (VD)  and  upper  mantle  (Vc) 


Fia,  2:  Schematic  diagram  of  wave  propagation 
through  the  crust-mantle  system 


Furthermore  it  had  been  found  necessary  to 
divide  the  lower  crust  into  several  gradient 
layers  in  order  to  account  for  the  more  im¬ 
portant  arrivals  seen  on  the  seismograms* 

From  a  comparison  of  seismic  refraction  and 
reflection  data  in  Central  Europe  it  was 
concluded  that  the  next  two  interfaces  just 
below  the  sialic  low-velocity  zone  seem  to 
be  associated  with  the  so-called  "Conrad 
discontinuity"  in  the  northern  and  southern 
parts  of  Germany ,  respectively* 

Some  observations  provide  strong  support  to 
the  proposed  existence  of  an  intermediate 
crustal  layer  just  above  the  "Mohorovi£i6 
discontinuity"  with  velocities  of  about 
7*5  km/sec  (ref.  LANDISMAN  and  MUELLER,  1966; 
FUCHS  and  LANDISMAN,  1966  a+b)  .  This  high- 
velocity  layer  has  been  found  now  in  many  parts 
of  Europe*  Its  average  thickness  must  be  of  the 
order  of  10  km.  A  typical  case  will  be  pre¬ 
sented  later  on  which  emphasizes  the  neces¬ 
sity  of  introducing  a  second  velocity  reversal 
right  on  top  of  this  high-speed  layer. 


1*2  Wave  Propagation  through  the  Crust-Mantle  System 

If  the  validity  of  a  linear  theory  qf  wave  pro¬ 
pagation  along  a  ray  path  from  the  source  through 
the  crust  and  mantle  to  the  receiving  station 
is  assumed,  the  wave  forms  of  the  P  and  S  signals 
at  large  enough  distances  are  shaped  by  a  number 
of  factors  which  depend  on  the  history  of  the 
signal  as  indicated  in  Fig*  2  According  to  this 
schematic  diagram  it  should  be  possible  to  trace 


back  the  signal  character  of  a  certain  phase 
in  the  seismogram  at  the  receiving  station 
through  the  recording  instrument*  the  receiver 
system*  the  "homogeneous  mantle"  and  the  source 
system  to  the  source. 

In  a  realistic  model  of  the  earth  the  complex 
structure  of  the  crust  and  upper  mantle  both 
at  the  source  and  the  receiver  location  must 
be  taken  into  account.  To  a  first  approximation 
the  inhomogeneous  crust  and  uppermost  part  of 
the  mantle  are  replaced  by  a  system  of  (n-1) 
homogeneous  isotropic  layers  with  plane  parallel 
interfaces  as  shown  in  Fig.  2  (ref.  FUCHS*  1966). 
The  compressional  velocity*  the  shear  velocity* 

A.U 

the  density  and  the  layer  thickness  of  the  i— 
layer  have  been  denoted  by  ot>  (■  Vpi)  * 

(-  Vgi)  *  Si*  and  ^i#  rc8Pectively*  Primed 
symbols  refer  to  the  "receiver  system"*  while 
unprimed  quantities  describe  the  "source  system". 
The  top  layer  (index  1)  in  either  of  the  two 
systems  is  specified  to  have  a  free  surface. 

The  layer  representing  the  deeper  mantle  is 
assumed  to  be  relatively  homogeneous*  i.e.  the 
elastic  properties  in  this  "half  space"  vary 
more  gradually  than  they  do  in  the  differentiated 
system  of  the  crust  and  upper  mantle.  In  Fig.  2 
the  angle  is  the  angle  of  incidence  for  P 
waves  at  the  base  of  the  receiver  system  which 
for  a  particular  ray  is  equal  to  the  angle  of 
radiation  into  the  mantle  at  the  base  of  the 
source  system. 

A  vertical  point  force  acting  on  the  free  surface 
of  the  source  system  has  been  used  as  the  simplest 
model  of  an  atmospheric  explosion.  For  underground 
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explosions  an  Impulsive  source  situated  at  a 
shallow  depth  h  In  the  top  layer  has  been 
assumed*  If  the  dimensions  of  the  source 
region  are  small  compared  to  the  distance  from 
the  layer  boundaries ,  the  source  may  be  regarded 
as  a  concentrated  point  source* 


2.  THE  RECEIVER  SYSTEM 

Before  the  distorting  effects  of  the  earth's 
crust  and  upper  mantle  on  the  spectrum  of 
seismic  signals  are  discussed,  It  should  be 
emphasized  that  all  signal  phases  recorded  on 
a  seismogram  have  to  be  corrected  for  the 
amplitude  and  phase  response  of  the  recording 
instrument*  In  the  subsequent  analysis  proper 
care  is  taken  of  this  instrumental  correction 
except  for  cases  where  it  is  specifically 
mentioned. 

The  "Receiver  Crust" 

The  comparatively  inhomogeneous  structure  of 
the  crust  and  upper  mantle  gives  rise  to 
considerable  signal  distortion  as  has  been 
pointed  out  before*  The  spectral  behavior  of 
the  receiver  system  must  therefore  be  known 
quite  accurately,  and  this  in  turn  requires  a 
detailed  knowledge  of  crustal  and  upper  mantle 
structure  in  the  European  area.  A  comprehensive 
survey  of  crustal  investigations  in  Central 
Europe  has  been  undertaken  in  order  to  provide 
the  required  information* 
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In  Fig.  3  all  seismic  refraction  profiles  are 
compiled  which  have  been  studied  since  1947, 

For  about  ten  years  all  geophysical  institutes, 
state  geological  surveys  and  geophysical  prospec¬ 
ting  companies  in  the  Federal  Republic  of  Germany 
cooperate  closely  within  an  extensive  research 
program  called  the  "Determination  of  Crustal 
Structure  in  Central  Europe".  This  program  is 
sponsored  by  the  German  Research  Association 
(Deutsche  Forschungsgemeinschaft) * 

Whorever  possible  observations  were  carried  out 
in  reverse  directions  and  along  overlapping 
profiles.  These  are  marked  in  Fig.  3  by  double 
lines  and  broken  lines,  respectively.  The 
shotpoints  are  denoted  by  black  dots  at  the 
end  of  each  line  and  are  numbered  consecutively. 

A  preliminary  interpretation  of  most  of  these 
measurements  has  been  published  some  years  ago 
by  the  GERMAN  RESEARCH  GROUP  FOR  EXPLOSION 
SEISMOLOGY  (1964). 

.2  Travel  Times  of  P  Waves  and  the  Velocity 
Structure  of  the  Crust 

(G.WOLBER) 

An  example  for  a  presumably  homogeneous  path 
of  propagation  is  provided  by  a  line  running 
southeastward  from  shotpoint  10  (VOGGENDORF) 
in  Fig.  3  parallel  to  the  margin  of  the 
Bohemian  Massif.  This  profile  starts  about 
90  km  east  of  the  Gr&fenberg  Seicmological 
Observatory  in  the  "Moldanubicum" ,  It  is 
240  km  long  and  extends  far  into  Austria.  A 
total  of  52  stations  with  identical  instruments 
were  occupied,  all  located  on  crystalline 
outcrops  of  the  Bohemian  Massif. 
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In  Fig.  4  a  seismogram  section  with  correlations 
is  shown  for  this  profile  VOGGENDORF  -  SE.  All 
records  were  digitized  manually  and  then  brought 
to  the  samo  reduced  time  scale  by  a  digital 
plotter.  Four  major  travel-time  branches  (a,  b, 
c,  d)  could  be  identified.  A  velocity-depth 
structure  Vp  ( z )  has  been  determined  from  the 
correlations  drawn. 

The  result  is  shown  in  Fig,  7  (left) .  Delow  a 
pseudo-gradient  zone  near  the  surface  two 
velocity  reversals  had  to  be  introduced  in 
order  to  explain  the  overlapping  travel-time 
curves  b  and  c  in  Fig.  4  (ref.  LAMDIS’lAN, 

MUELLER  and  FUCHS,  1967) .  As  mentioned  earlier 
the  existence  of  a  high-speed  layer  with  a  P 
velocity  of  7.8  km/sec  just  above  the  classical 
MohoroviciC  discontinuity  must  be  assumed  if  the 
duplication  of  the  first  arrivals  at  greater 
distances  is  to  be  explained.  These  new  crustal 
features  are  presently  explored  in  more  detail. 

Amplitude  Decay  and  the  Qp  Structure  of  the  Crust 

(C.WOLBER) 

Following  the  schematic  diagram  in  Fig,  2  the 
data  processing  starts  with  a  spectral  analysis 
of  the  observed  seismogram  which  is  then  corrected 
for  the  effects  of  instrumental  distortion.  A 
schematic  flow  diagram  of  the  harmonic  analysis 
procedure  is  presented  in  Fig,  6.  The  computer 
program  reads,  checks,  filters,  decimatos  and 
windows  the  data  before  it  performs  the  frequency 
analysis.  Optionally  the  output  is  listed,  plotted 
or  punched  on  cards. 


INPUT: 

1.) 

DATA  CONTROL  CARD 

2.) 

DATA  CARDS 

3.) 

CASE  PARAMETER  CARD 

4.) 

END  OF  CASE  CARD 

OPTIONALLY: 

5.) 

FILTER  WEIGTIiS 

6.) 

TIME  WINDOW  WEIGTHS 

SAFETY 

PACKAGE: 

ciie 

CKING  FOR* 

l.) 

ORDER  OF  CARDS 

2.) 

NUMBER  OF  DATA  POINTS 

3.) 

PUNCHING  ERRORS 

4.) 

END  OF  CASE  CARD 

DATA  PREPARATION 


1. )  SHIFTING  OF  DATA 

2. )  REMOVAL  OF  LINEAR  TREND 

3. )  REMOVAL  OF  AVERAGE 
OUTPUT  OF  TREND  AND  AVERAGE 


FREQUENCY  ANALYSIS 


OUTPUT* 


1. )  LISTING  OF 

FREQUENCY,  AMPLITUDE,  PHASE 

2. )  PLOTTING  OF 

AMPLITUDE,  PHASE 

3. )  OPTIONAL  OUTPUT  ON  CARDS 


6*  Schematic  flow  diagram  of  spectral  analysis  program 


Spectral  analyses  of  the  correlated  seismogram 
phases  in  the  record  section  of  Fig,  4  have 
been  used  to  deduce  an  average  specific  quality 
factor  Q  for  P  waves  (Qp)  in  the  crust  of  the 
Bohemian  Massif.  If  the  amplitude  response 
of  the  explosive  source  and  of  the  receiver 
are  known#  any  systematic  variation  of  the 
observed  signal  spectrum  must  be  due  to  variation 
in  the  spectral  response  of  the  propagation 
system.  This  response  is  determined  by  the 
geometry  of  the  transmission  path  and  by 
anelastic  deficiencies  of  the  propagation  medium. 
For  distances  large  compared  with  the  depth  of 
penetration  the  geometry  should  not  have  any 
influence  on  the  amplitudes  unless  the  medium  is 
severely  inhomogeneous  in  lateral  directions.  Under 
these  simplifying  assumptions  any  frequency- 
dependent  amplitude  decay  may  be  attributed  to 
absorption  and  -  to  a  minor  degree  -  to  scattering. 

In  Fig.  5  the  amplitude  spectra  are  shown  for  the 
travel-time  branch  c  (7.8  km/soc)  as  correlated 
in  Fig.  4.  All  the  spectra  in  our  analysis  have 
been  corrected  individually  for  instrumental 
distortion.  Following  a  well-known  procedure 
(see  e.g.  O'BRIEN#  1967)  the  spectral  gradients 
were  then  determined  for  each  spectrum  in  the 
distance  range  given.  From  a  plot  of  spectral 
gradients  versus  range  an  average  value  for  Qp 
could  be  estimated  provided  a  representative 
value  for  the  P  wave  velocity  (Vp)  was  chosen. 

Assuming  ideal  head-wave  propagation  the  velo¬ 
city-depth  distribution  in  Fig,  7  (left)  could 
be  used.  With  this  model  a  least-squares  solution 
led  to  a  mean  Qp  of  670  —  100  for  the  upper 


crust  and  of  670  -  150  for  the  lower  crust 
(see  Fig.  7 ,  right)  .  For  the  uppermost  mantle 
a  Q_  of  250  -  80  has  been  found.  There  are 
some  indications  that  a  gradual  decrease 
of  Qp  occurs  in  the  transition  region  between 
crust  and  mantle.  Determinations  of  Qp  for 
the  middle  crust  gave  values  around  450  with 
a  relatively  high  error.  The  Qp  structure  presen¬ 
ted  in  Fig.  7  (right)  is  far  from  being  complete, 
but  demonstrates  that  the  Qp  values  for  the 
crust  inferred  from  Qg  measurements  (ref.  PRESS, 
1964;  ANDERSON,  BEN-MENAIIEM  and  ARCHAMBEAU,  1965) 
seem  to  be  too  high  by  a  factor  of  about  two. 

Crustal  Transfer  Functions 

(K.  FUCHS) 

Since  results  of  seismic  refraction  measurements 
are  only  available  for  a  few  areas  in  Europe 
a  program  has  been  initiated  to  check  and 
supplement  these  data  by  determining  the 
influence  of  the  "receiver  crust"  on  the 
spectral  behavior  of  long-period  body  waves. 
Theoretical  transfer  functions  for  the  crust- 
mantle  system  have  been  calculated  following 
earlier  investigations  by  PHINNEY  (1964)  , 

HANNON  (1964),  FERNANDEZ  (1967)  and  BEN-MENAIIEM, 
SMITH  and  TENG  (1965)  •  A  comparison  of  these 
theoretical  transfer  functions  with  experimen¬ 
tal  results  have  enabled  us  to  test  models  of 
crustal  structure  in  regions  whore  the  results 
of  other  methods,  such  as  travel  time  studies 
and  surface-wave  dispersion  measurements,  did 
not  lead  to  a  consistent  picture. 
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Fig,  8 1  Comparison  of  experimental  and  theoretical 

crustal  transfer  functions  for  Denmark  (C0n) . 
i  is  the  angle  of  incidence  at  the  base  of 
the  crust. 
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Depth  Thickness  Vp  Vg  Density 
(km)  (km)  (km/s)  (km/s)  (g/cm3) 


1.9 

8.3 

27.8 

1.9 

6.4 

19,5 

Mode  1 

DENS 

20  -  Denmark 

1.9 

8.3 

29.3 

1.9 

6.4 

21.0 

Model 

DENS 

21  -  Denmark 

2.70 

11.23 

29.40 

2.70 

8.53 

18.17 

Model 

JTLN 

01  -  Denmark 

5.10 

11.77 

30.00 

5.10 

6.67 

18.23 

Model 

JTLM 

01  -  Denmark 

1.05 

12.20 

32.80 

1.05 

11.15 

20.60 

Model 

JTLS 

01  -  Denmark 

3.0  1.732  2.33 
6.1  3.522  2.86 
6.6  3.811  3.00 
7.9  4.561  3.38 


(HIRSCHLEEER,  HJELME, 
SELLEVOLL,  1966) 


3.0 

1.732 

2.33 

6.1 

3.522 

2.86 

6.6 

3.811 

3.00 

8.2 

4.734 

3.43 

(IIIRSCHLEBER, 

IIJELME, 

SELLEVOLL,  1966) 

3.0 

1.73 

2.20 

6.0 

3.46 

2.75 

6.7 

3.87 

2.95 

8.1 

4.67 

3.40 

(IIIRSCHLEBER, 

1964; 

pers 

.  communication) 

3.0 

1.73 

2.20 

6.0 

3.46 

2.75 

6.7 

3.87 

2.95 

8.1 

4.67 

3.40 

(IIIRSCHLEBER, 

1964; 

pers 

,  communication) 

3.0 

1.73 

2.20 

6.0 

3.46 

2.75 

6.7 

3.87 

2.95 

8.1 

4.67 

3.40 

(IIIRSCHLEBER,  1964; 
pers,  communication) 


Table  It  Crustal  Models  in  Denmark 


Depth 

Thickness 

VP 

vs 

Density 

(km) 

(km) 

(km/s) 

(km/s) 

3 

(g/cm 

4.75 

21.75 

33.00 

4.75 

3.98 

2.38 

2.34 

17.00 

5.90 

3.30 

2.82 

11.25 

6.58 

3.80 

2.92 

8.14 

4.70 

3.30 

Model 

IB 

1 

-  Iberian 

Peninsula  (PAYO, 

1965) 

A  7 

4.7 

3.94 

2.37 

2.34 

23.0 

33.0 

18.3 

5.90 

3.30 

2.80 

10.0 

6.60 

3,70 

2.85 

8.10 

4.72 

3.30 

Model 

IB 

2 

-  Iberian 

Peninsula  (PAYO, 

1965) 

Table  2t 


Crustal  Models  of  the  Iberian  Peninsula 
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As  an  example  the  spectral  ratio  TE  (y) 
of  the  vertical  component  W  (V)  to  the 
horizontal  component  U  (V)  is  shown  in  Fig*  8 
f jr  two  deep-focu3  earthquakes  as  recorded 
at  Copenhagen  (COP) ,  Denmark.  The  theoretical 
transfer  functions  TEN  ( V)  for  various  crustal 
models  of  Denmark  are  compared  with  the 
experimental  results*  It  is  seen  that  the 
crustal  models  DENS  20  and  JTLN  01  (see  Table  1) 
come  closest  to  the  observations  with  regard  to 
peak  ratio  and  frequency. 

A  similar  analysis  has  been  carried  out  for 
station  Toledo  (TOL) ,  Spain.  It  could  be  shown 
that  the  crustal  model  ID  2  (see  Table  2)  of 
PAYO  (1965)  provided  a  much  better  fit  to  the 
data  than  model  IB  1  (ref.  Table  2)  which  we 
had  used  in  our  previous  investigations 
(MUELLER  and  FUCHS,  1966). 

3.  THE  SOURCE  SYSTEM 

(K.  FUCHS) 

3.1  The  "Source  Crust" 

Since  no  direct  determinations  of  crustal 
structure  are  available  for  the  source  areas 
covered  in  our  studies  the  elastic  properties 
of  the  earth's  crust  in  those  regions  had  to 
be  inferred  by  extrapolation  of  published  re¬ 
sults.  For  Novaya  Zemlya  and  the  Barents  Sea 
the  two  crustal  models  NUR-1  (PENTTIL&,  1965) 
and  URAL  (DEMENITSKAYA,  1961)  depicted  in 
Fig.  9  were  chosen  to  give  an  estimate  of  the 
crustal  structure  under  Novaya  Zemlya* 


NURMIJ  ARVI  NUR-1 
(PENTTILAJ965) 


URAL 

(DEMEN  I TSKAYA,  1961) 

‘i 

i 

i 

i 


*• 1 


9:  "Source  crust"  models  NUR-1  and  URAL 


»  (mF]j  [TT]*  EZ3*  Elil* 


10:  Crustal  section  for  the  profile  UCHTA-KEM 
in  Karelia  {U,S,,S,R,)  -  Explanation  of 
legend  see  text 
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A  better  approximation  is  probably  given  by  the 
crustal  section  of  LITVINENKO  and  NEKIIASOVA 
(1962)  for  the  profile  UCHTA-KEM  in  Karelia 
reproduced  in  Fig.  10.  The  symbols  used  in 
the  section  are  explained  as  follows t 

1  “  Interface  found  from  head  waves 

2  =>  Interface  found  from  reflected  waves 

3  «=  Fault  zone  determined  by  deep  seismic  sounding 

4  =  Fault  zone  determined  by  field  geology 

5  «»  Granite  bodies 

6  =  Granodiorite  zones 

7  =  Regions  with  basic  and  ultrabasic  rocks 

The  values  of  P  velocity  are  given  in  meters  per 
second.  It  is  seen  that  the  main  features  of  this 
crustal  model  are  similar  to  the  ones  of  our 
preferred  model  URAL.  Both  structures  are 
presently  investigated  if  they  do  provide  a 
satisfactory  explanation  for  the  spectral  behavior 
of  long-period  body  waves  recorded  at  nearby 
observatories  in  Finland. 

Source  Mechanisms 

A  vertical  point  force  acting  on  the  free  surface 
of  a  layered  system  (see  Fig.  2)  ha3  been  used  in 
an  attempt  to  evaluate  the  influence  of  the  earth's 
crust  and  the  uppermost  mantle  on  the  spectrum 
of  seismic  body  waves  excited  by  atmospheric 
explosions.  The  radiation  pattern  for  the  far 
field  has  been  calculated  at  a  number  of  discrete 
frequencies  for  the  two  crustal  models  NUR-i  and 
URAL  (see  Fig.  9).  Fig,  11  demonstrates  how 
differences  in  crustal  structure  of  the  source 
region  affect  the  shape  of  the  radiation  pattern. 


V- 
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Fig,  11:  Radiation  pattern  for  a  vertical  point  source  acting 
on  the  free  surface 
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For  underground  explosions  an  impulsive  point 
source  buried  at  some  depth  within  the  crust 
served  as  an  adequate  model.  The  transfer 
function  for  thi3  configuration  has  been 
derived  by  FUCHS  (1966).  In  Fig.  12  the  ra¬ 
diation  patterns  in  the  frequency  range  from 
0  to  0.5  cps  are  shown  for  the  two  "source 
crust"  models.  As  before  only  the  right  half 
of  the  symmetrical  pattern  has  been  drawn.  A 
clear  directional  effect  is  seen. 

The  radiation  pattern  of  the  buried  explosive 
source  differs  significantly  from  that  of  the 
vertical  point  force.  In  the  case  of  the  buried 
source  a  strong  primary  reflection  is  generated 
at  the  free  surface  of  the  crust.  Its  inter¬ 
ference  with  the  primary  wave  strongly  shapes 
the  radiation  pattern.  The  surface  reflection 
is,  of  course,  absent  in  the  case  of  the  point 
force  acting  on  the  free  surface.  Experimentally 
determined  amplitude-distance  curves  will  be 
influenced  by  the  frequency-dependent 
distortions  of  the  radiation  pattern.  This  effect 
is  presently  being  analyzed  more  thoroughly. 
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Fig,  12:  Radiation  pattern  for  an  explosive  point  source 
buried  at  a  depth  cf  0,5  km 
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4.  STRUCTURE  OF  THE  UPPER  MANTLE 
IN  EUROPE 


4.1  Propagation  of  Seismic  Body  Waves 

4.11  Travel  Times  from  Large  Explosions 

(J.  ANSORGE) 

Travel  time  studies  of  earthquake  events  suffer 
from  the  inherent  inaccuracies  of  epicenter 
location  and  origin  time.  We  have  therefore 
tried  to  obtain  more  accurate  travel  time  data 
in  Europe  for  the  distance  range  to  10°  by  using 
large  explosions.  A  number  of  artificial  events 
recorded  along  suitable  long-range  profiles  were 
carefully  analyzed  in  an  attempt  to  resolve  the 
apparent  discrepancies  between  earthquake  and 
explosion  observations.  In  Fig.  13  long-range 
seismic  refraction  lines  are  shown  for  4  selected 
events i 

(1)  The  HELGOLAND  Explosion  of  18  April  1947 

Charge  size  about  8000  tons  of  ammunition 
Length  of  SSE  profile!  640  km 

(2)  The  LAC  NEGRE  Explosions  of  September  1966 

Charge  sizes  up  to  40  tons 
Length  of  N  profile!  700  km 

(3)  The  LAC  de  l'EYCHAUDA  Explosion  of  22  September 

1963  (BEAUFILS,  MECHLER  and  ROCARD,  1965) 
Charge  size!  2  tons 

These  data  could  be  combined  with  obser¬ 
vations  of  (2)  to  provide  an  overlapping 
profile. 

Length  of  MW  profile!  800  km 

(4)  The  FOLKSTONE  Explosion  of  22  July  1967 

Charge  size  about  500  tons 
(Magnitude  M  »  4.7  according  to  USCGS) 
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F  -  FOLKSTONE  July  22,1967 
H  =  HELGOLAND  April  18, 1947 
LE-LAC  DE  L’EYCHAUDA  Sept.  1963 
LN-LAC  NEGRE  Sept.1966 


oiauons  wmch  recorded 
FOLKSTONE  Explosion  in  two  main 
directions 


Refraction  Profiles 


Map  of  long-range  seicr.ic  refraction  profiles 
in  Europe 
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(a)  Observations  in  Central  Europe  to  10°  - 
The  earthquake  stations  in  Central  Europe 
were  groped  in  two  main  azimuthal 
directions  -  a  distinction  which  turned 

out  to  be  unnecessary  in  this  distance  range 

(b)  Observations  in  Scandinavia  to  20°  - 
presently  being  analyzed. 

All  these  long-range  seismic  measurements  have 
shown  that  the  distribution  of  compressional 
velocity  (Vp)  with  depth  in  the  uppermost  part 
of  the  mantle  is  much  more  complicated  than 
hitherto  assumed.  It  also  differs  significantly 
from  region  to  region.  For  a  long  time  it  has 
been  suspected  that  a  velocity  inversion  for 
P  waves  should  accompany  the  shear-wave  low- 
velocity  channel  found  in  the  upper  mantle  at 
depths  between  about  100  and  200  km. 

The  timed  explosions  so  far  analyzed  have  con¬ 
sistently  produced  evidence  for  a  triplication 
of  the  P  wave  travel  time  curve  in  the  distance 
range  between  4°  and  about  7°.  With  the  high 
time  resolution  of  our  new  field  equipment  a 
new  phase  could  be  identified  which  was  termed 
Pm,  (see  Fig.  14) ,  Typical  values  for  the  sur¬ 
face  velocity  of  this  phase  lie  between  8.25 
and  8,35  km/sec,  while  for  Pn  (=  P^)  relatively 
low  values  between  8,05  and  9,10  km/sec  are 
measured.  The  new  phase  explains  the  previously 
observed  break  or  jump  in  the  first  arrival 
curve  around  9°. 

In  Fig.  15  the  explosion  data  of  the  LN-LE  MW 
profile  (see  Fig,  14)  are  plotted  in  conjunction 
with  earthquake  data  for  events  of  about  the 
same  epicentral  location.  If  the  latter  are 
reduced  to  zero  focal  depth  both  sets  of  ob¬ 
servations  can  be  combined  yielding  a  very 


s, 


Explosions 

Earthquakes 


Fig,  1 5 1  Travel  times  from  large  explosions  and  oarthnu  i):cg  in 
Central  Europe  (P  -  Pd)  with  deduced  velocity-depth 
model  13 
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consistent  travel  time  pattern  for  Central  Euro¬ 
pe. 

The  velocity  distribution  for  Model  ©  depicted 
in  Fig.  15  provides  a  satisfactory  explanation 
of  the  data.  Two  velocity  reversals  at  depths 
between  70  and  220  km  are  required  in  the  model 
in  order  to  account  for  the  observed  overlap 
of  the  various  travel-time  branches  (ref. 
MUELLER,  ANSORGE,  MAYER-ROSA  and  FUCHS,  1968). 

Additional  information  on  the  structure  of  the 

uppermost  mantle  in  Central  Europe  was  obtained 

from  the  recordings  of  the  FOLKSTONE  Explosion 

(22  July  1967)  .  The  location  of  this  event  is 

precisely  known,  but  the  origin  time  given  by 

USCGS  and  BCIS  differ  by  1.3  seconds.  This 

discrepancy  could  be  resolved  by  comparing  the 

P  observations  with  recent  seismic  refraction 
n 

measurements  in  the  British  Channel  between 
Cornwall  and  the  Bretagne.  From  the  apparent 
velocity  and  intercept-time  the  USCGS  deter¬ 
mination  was  found  to  be  correct. 

In  Fig,  16  all  readable  signal  onsets  are  plot¬ 
ted  as  a  function  of  distance.  Three  travel¬ 
time  branches  could  be  correlated  which  seemed 
to  be  independent  of  azimuth.  The  duplication 
of  the  first  arrivals  in  the  distance  range 
between  4°  and  5°  led  to  the  identification 

of  a  P  .  and  a  P-  branch  with  velocities  of 
ml  m2 

8,12  and  8.31  km/sec,  respectively.  At  close- 

in  distances  a  P^  velocity  of  8.05  km/sec  was 

measured.  The  termination  of  the  P^  line 
Q  n 

around  5  for  all  profiles  in  Europe  and  the 

relatively  short  critical  distances  for  the 

P„,  and  P  -  branches  has  forced  us  to  intro- 
mi  m2 

duce  an  additional  low-velocity  zona  between 
45  and  60  km  depth. 


-IOA(scc) 
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These  results  confirm  the  complexity  of  the 
uppermost  mantle  and  indicate  that  at  least 
three  velocity  reversals  occur  between  the 
Mohorovicic  discontinuity  and  the  abrupt 
velocity  increase  postulated  by  LEHMANN 
(1959)  at  a  depth  of  about  220  km,  namely 


between 

45 

and 

60 

km, 

80 

and 

100 

km, 

140 

and 

220 

km. 

4*12  Travel  Times  and  Amplitude  Snectra 

from  Earthquakes  (D.  MAYER-ROSA) 

At  larger  distances,  i.e.  for  5°<  30°, 

the  travel-time  data  were  supplemented  by 
earthquake  observations  which,  of  course, 
lack  the  time  resolution  of  the  explosion 
measurements.  In  addition  they  are  subject 
to  the  usual  uncertainties  in  the  epicenter 
and  origin  time  determinations.  A  total  of 
600  events  recorded  at  Stuttgart  (STU)  bet¬ 
ween  1937  and  1967  was  chosen  for  this  3tudy, 
Epicentral  distances  were  computed  based  on 
the  focal  data  of  the  I.SS  up  to  1961,  and 
using  the  USCGS-Bulletin  for  the  period  1961 
through  1967.  Accordingly  the  standard  (root- 
mean-square)  error  of  the  travel-time  deter¬ 
minations  for  P  waves  is  +  1,5  sec  between 
1937  and  1960,  and  +  0,5  sec  between  1961 
and  1967.  For  S  waves  the  error  is  about 
twice  as  high. 

In  Figs,  17  and  18  the  measured  P  and  S 
travel  times  are  shown  for  two  profiles 
directed  SW  and  SE  of  Stuttgart,  i.  e.  towards 
Spain  and  Greece.  They  were  determined  by 
observations  of  a  multitude  of  events  at  a 
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17 i  P  wave  travel  times  from  earthquakes  (h=0) 


—  SW 

—  SE 


10  15  20 
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18 i  S  wave  travel  times  from  earthquakes  (h“0) 
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fixed  station  (e.g.  STU)  and  of  selected  events 
of  suitable  magnitude  along  a  chosen  array  of 
European  stations.  An  additional  check  was  nro- 
vided  by  the  observed  travel  times  of  inter¬ 
mediate  and  deep  focus  earthquakes.  All  travel¬ 
time  observations  are  corrected  to  zero  focal 
depth  and  are  reduced  by  10A  for  P  and  20 A 
for  S  waves.  In  both  cases  five  travel-time 
branches  can  be  seen,  namely 


■  p 

S,  -  S 

n 

d 

P 

S 

m 

m 

P 

s 

rl 

rl 

Pr2 

Sr2 

Ptele 

Stej.e 

The  apparent  velocities  of  the  first  two  branches 
are  different  for  the  SW  and  SE  profiles,  while 
the  other  three  do  not  show  any  measurable 
differences  in  velocity.  The  SW  travel-time 
curves,  however,  are  consistently  delayed  by 
2  sec  in  P,  and  by  about  3  sec  in  S  compared  to 
the  SE  line.  This  immediately  demonstrates  that 
there  must  be  pronounced  regional  differences 
in  upper  mantle  structure.  The  considerable 
overlap  of  the  various  travel-time  branches 
indicates  that  a  much  higher  velocity  contrast 
must  exist  at  the  corresponding  depth,  much 
higher  than  is  suggested  by  the  measured  apparent 
velocities  at  the  surface.  Distinct  velocity  re¬ 
versals  for  some  depth  intervals  are  therefore 
required  to  reconcile  the  results  of  theoretical 
computations  with  the  observations. 


In  Fig.  19  examples  of  the  P  and  S  wave  groups 
are  shown  as  recorded  at  distances  of  8?2, 

13?6  and  16?7  (SE  of  Stuttgart) »  The  magnitudes 
range  from  6.0  to  6.5.  As  the  recording  instru¬ 
ment  (Galitzin  12-12)  and  the  "receiver  crust" 
is  the  same  in  all  cases,  and  if  it  is  assumed 
that  the  "source  crust"  is  about  the  same  and 
that  there  is  no  directional  effect  of  the 
source  for  the  periods  observed,  then  the 
spectral  content  of  the  various  phases  of  the 
signal  is  solely  determined  by  the  propagation 
medium,  i.e.  by  the  path  through  the  rather 
complex  structure  of  the  upper  mantle. 

The  dominating  periods  (T)  of  P  and  S  wave 
arrivals  were  determined  as  a  function  of 
distance  by  Fourier  analysis  and  by  visual 
analysis  of  the  seismogram  phases,  as  well  as 
by  spectrograph  techniques  (ref.  EWING,  MUELLER, 
LANDISMAN  and  SATO,  1959) »  The  results  are  shown 
in  Figs.  20  and  21.  A  more  or  less  pronounced 
absorption  effect  can  be  seen  for  each  of  the 
different  phases.  Average  Q  values  could  be 
estimated  from  these  data  which  are  difficult 
to  interpret  in  view  of  the  complex  velocity 
structure.  The  dominating  periods  increase 
systematically  as  the  energy  penetrates  deeper 
and  deeper  into  the  mantle.  Except  for  the 
♦•elesoismic  P  and  S  branches  the  predominant 
periods  in  S  are  nearly  twice  as  long  as  the 
values  for  P.  From  these  periods  and  the 
corresponding  wave  lengths  a  lower  limit  for 
the  thickness  of  the  various  layers  through 
which  the  waves  travel  could  be  deduced. 


o 


DOMINATING  PERIODS 
OF  P  ARRIVALS 
GALUZIN  12-12 


o  Fourier  Analysis 
1  Visual  Analysis 


Dominating  wave  periods  of  P  arrivals 
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21 i  Dominating  wave  periods  of  S  arrivals 
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By  trial  and  error  velocity-depth  distributions 
were  derived  which  can  explain  the  observed  P 
and  S  travel  times  given  in  Figs.  17  and  18, 

They  are  depicted  in  Fig,  22*  From  the  argu¬ 
ments  presented  above  and  similar  conclusions 
arrived  at  by  DOWLING  and  NUTTLI  (1964)  it  was 
necessary  to  introduce  a  number  of  velocity 
reversals  in  order  to  explain  the  overlapping 
travel-time  segments.  Pronounced  regional 
differences  between  the  SW  and  SE  profiles  are 
noted  down  to  depths  of  about  320  km*  It  should 
be  pointed  out  specifically  that  a  third  low- 
velocity  channel  in  S  with  a  lower  velocity 
contrast  had  to  be  assumed  in  the  model  calcu¬ 
lations.  The  Srl  -  Sf2  travel-time  branches 
can  only  be  explained  by  a  velocity  inversion 
between  260  and  320  km  depth.  For  P  waves  the 
corresponding  channel  is  either  absent  or  so 
weak  that  it  cannot  be  resolved  with  the 
present  means  of  observation. 

In  Fig,  22  the  velocity-depth  profiles  for 
Europe  as  proposed  by  LEHMANN  (1959,  1961) 
are  given  for  reference  and  comparison.  It  can 
be  seen  that  the  new  P  distribution  oscillates 
around  the  old  average  value  down  to  a  depth 
of  about  210-220  km,  A  similar  "lamellar" 
structure  of  the  upper  mantle  is  found  for 
S.  The  velocity  inversions  are  much  more  severe 
than  the  single  reversal  which  has  been  assumed 
to  exist  so  far.  If  the  lower  channel  boundaries 
are  indeed  as  sharp  as  LEHMANN  (1959,  1961) 
and  others  have  suggested,  then  clear 
reflections  should  be  observed.  Evidence  is 
found  in  the  literature  for  North  America 
(ref.  KNOPOFF,  MUELLER  and  PILANT,  1966)  which 
confirms  the  postulate  of  LEHMANN  (1962,  1964) • 
In  Europe  we  were  able  to  trace  at  least  parts 
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of  two  reflection  hyperbolas  which  can  be 
associated  with  pronounced  discontinuities  at 
depths  of  about  210  and  540  km. 

Supporting  evidence  for  the  velocity-depth 
profiles  presented  in  Fig,  22  comes  from  the 
distribution  of  earthquake  focal  depths  in 
Europe.  In  Fig,  23  the  number  of  earthquakes  is 
plotted  as  a  function  of  focal  depth  (h)  for  the 
time  interval  1960-1967.  Theso  data  are  more 
reliable  than  the  ones  published  previously  for 
the  European  area.  If  the  average  depth  ranges 
of  the  zones  of  velocity  reversals  are  compared 
to  the  focal  depth  distribution  it  is  found  that 
they  correlate  well  with  the  observed  relative 
minima.  In  particular  the  concentration  of  foci 
close  to  the  upper  boundaries  of  the  low- 
velocity  channels  as  well  as  the  complete 
absence  of  earthquakes  with  focal  depths  around 
200  km  should  be  noted.  It  is  thus  strongly 
emphasized  that  the  structure  of  the  upper 
mantle  is  much  more  complex  than  was  implied 
by  our  previous  models. 

4 . 2  Dispersion  of  Seismic  Surface  Waves 

(D.  SEIDL) 

Extensive  studies  of  the  dispersion  character¬ 
istics  of  longer  period  surface  wavos  of  the 
Love  and  Rayleigh  type  have  confirmed  the 
existence  of  low-velocity  layers  in  the  upper 
mantle.  They  have  also  shown  that  there  are 
pronounced  regional  variations  in  the  velocity 
distribution  with  depth.  The  attenuation  of 
long-period  surface  waves  and  the  decay  of 
the  free  oscillations  of  the  oarth  depend  in 
a  very  complicated  way  on  the  depth  distribution 
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of  parameters  of  elasticity  and  anelasticity . 
Complications  caused  by  the  source  mechanism, 
instrumental  distortion,  path  differences  as 
well  as  low-velocity  zones  are,  however,  con¬ 
siderably  reduced  compared  with  body  waves. 

In  Fig,  24  a  block  diagram  of  the  data-pro- 
ccssing  procedure  used  in  our  surface  wave 
studies  is  shown.  Starting  with  a  seismogram 
which  can  be  thought  of  being  a  superposition 
of  m  normal  (and  leaking)  modes  the  group  and 
phase  velocity  dispersion  as  v/oll  as  the  ab¬ 
sorption  coefficient  as  a  function  of  frequency 
(  V)  can  be  determined  for  any  two-3tation 
combination.  Presently  the  data  analysis  is  in 
progress  for  a  trans-European  profile  from 
KEVO  (Finland)  to  MALAGA  (Spain) . 

During  the  past  fev;  years  our  attention  has 
been  focussed  to  the  central  European  region 
and  most  recently  to  the  Rhinegrabcn  Rift- 
System  (MUELLER,  FETERSCHMITT,  FUCHS  and 
ANSORGE,  1967),  The  map  in  Fig,  25  summarizes 
the  surface-wave  observations  in  that  area. 
Heavy  lines  indicate  thb  station  network  used 
by  SEIDL  (1965)  and  SEIDL,  MUELLER  and  KNOPOFF 
(1966)  in  their  phase  velocity  measurements  to 
the  north  of  the  Alps.  Also  shown  by  solid 
lines  is  the  Alpine  network  of  KNOPOFF,  MUELLER 
and  PILANT  (1966),  Tripartite  results  in  the 
Alps  produced  evidence  for  a  well-developed 
low-velocity  channel  for  S  with  a  shear  velo¬ 
city  in  the  channel  of  4,2  km/soc.  Broken  lines 
in  Fig,  25  represent  the  new  network  around 
the  Rhinegraben  Rift-  System  which  is  nearly 
complete , 
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In  our  Annual  Summary  Report  No.  1  (1966)  a 
model  of  the  crust  and  upper  mantle  has  boon 
reviewed  which  was  termed  STUTTGART  ©  (rof . 
SCHNEIDER/  MUELLER  and  KNOPOFF,  1966).  It  was 
primarily  based  on  phase  and  group  velocity 
measurements  of  Rayleigh  waves  for  paths  in 
the  region  of  the  Alps  and  through  the  fore¬ 
land  to  the  north.  The  rather  scarce  body-wave 
travel  times  available  then  were  concordant  with 
that  model. 

Based  on  additional  studies  of  a  number  of  near 
earthquakes  (ref.  SCHNEIDER/  1964)  we  have  now 
deduced  a  new  model  which  seems  to  be  more 
representative  for  the  average  crustal  and 
upper  mantle  structure  in  Central  Europe.  This 
new  model/  called  STU-3  (see  Table  3  ) ,  consists 
of  a  four-layer  crust  in  agreement  with  seismic 
refraction  and  near-earthquake  observations.  No 
efforts  have  been  made  so  far  to  include  the 
two  low-velocity  zones  discussed  above  in  this 
average  crustal  model.  As  before  a  sharp 
discontinuity  is  postulated  at  a  depth  of  about 
220  km  which  forms  the  lower  boundary  of  the 
main  asthenosphere  low-velocity  channel. 

Observational  data  of  fundamental  mode  Rayleigh 
waves  for  the  lines  Stuttgart-Strasbourg  (SF.IDL/ 
MUELLER  and  KNOPOFF,  1966)  and  Stuttgart- 
Besan^on  (KNOPOFF,  MUELLER  and  PILANT,  1966) 
which  both  traverse  the  Rhinegraben  Rift  System, 
are  satisfactorily  explained  by  model  STU-3 
(see  Fig.  26) .  The  phase-velocity  results  for  the 
two  Alpine  lines  Stuttgart-Oropa  and  Stuttgart- 
Chur  are  in  better  agreement  with  the  old  model 
STUTTGART  1  (ref.  SCHNEIDER,  MUELLER  and 
KNOPOFF,  1966). 
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Depth 

(km) 


1.50 

4.50 
20.00 
30.00 

100.00 

220.00 

320.00 

410.00 

500.00 

600.00 

700.00 

800.00 

900.00 

1000.00 

1100.00 

1200.00 

1300.00 

1400.00 

1500.00 

1600.00 

1700.00 

1800.00 

1900.00 

2000.00 

2100.00 

2200.00 

2300.00 

2400.00 

2500.00 

2600.00 

2700.00 

2800.00 


Thickness  V  Vc  Density 


(km) 

(km/s) 

(km/s) 

(g/cm 

1.50 

3.50 

2.00 

2.50 

3.00 

5,60 

3.30 

2.75 

15.50 

6.00 

3.50 

2.85 

10.00 

6.70 

4.00 

3.00 

70.00 

8.15 

4.60 

3.30 

120.00 

8.20 

4.20 

3.40 

100.00 

8.49 

4.77 

3.53 

90.00 

8.81 

4.89 

3.60 

90.00 

9.32 

5.19 

3.76 

100.00 

9,97 

5.49 

4.01 

100.00 

10.48 

5.79 

4.23 

100.00 

10.85 

6.03 

4.41 

100.00 

11.12 

6.20 

4.55 

100.00 

11.33 

6.32 

4.64 

100.00 

11.49 

6,40 

4.71 

100.00 

11.64 

6.47 

4.77 

100.00 

11.78 

6  .53 

4.83 

100.00 

11.92 

6.59 

4.88 

100.00 

12.06 

6.65 

4.94 

100.00 

12.19 

6.70 

5.00 

100.00 

12.33 

6.76 

5.06 

100.00 

12.46 

6.81 

5.11 

100.00 

12.59 

6.85 

5.16 

100.00 

12.72 

6.90 

5.21 

100,00 

12.85 

6.95 

5.27 

100.00 

12.97 

7.00 

5.32 

100.00 

13.09 

7,05 

5.37 

100. GO 

13.21 

7  i  10 

5.42 

100.00 

13.33 

7.14 

5.47 

100.00 

13.46 

7.1b 

5.52 

100.00 

13.53 

7.23 

5.56 

100.00 

13.61 

7.28 

5.61 

13.64 

7.30 

5.66 

Table_3j  Crust  -  Mantle  Model  STU  -  3 
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Fie,  26:  Phase  velocity  measurements  of  Rayleigh  v;avcs  in 
entral  Europe  and  velocity-depth  distribution 
or  model  STU-3 
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In  contrast  to  this  previously  preferred 
structure  the  new  model  STU-3  contains  only 
one  layer  above  the  channel  in  the  upper  mantle, 
The  thickness  of  this  lid  (H^)  has  been  in¬ 
creased  by  20  km  to  70  km,  while  the  thickness 
of  the  channel  (»CH)  has  been  decreased  by  the 
same  amount  to  120  km.  All  shear  velocities 
above  the  channel  are  now  slightly  lower  than 
before. 

In  Fig,  27  the  partial  derivatives  of  phase  ve¬ 
locity  are  presented  for  model  STU-3,  The  upper 
half  of  Fig,  27  shows  how  the  phase  velocity 
(c)  of  Rayleigh  waves  is  altered  by  small 
variations  of  shear  velocity  in  the  channel 

,  In  the  lower  half  of  Fig,  27  the  changes 
of  c  caused  by  variations  of  the  lid  (H__)  and 

iiU 

channel  (HCJj)  thicknesses  are  displayed.  It  is 
seen  that  changes  in  the  channel  velocity  (flCH) 
are  an  order  of  magnitude  more  severe  than 
variations  in  the  geometrical  configuration ►  If 
for  the  time  being  the  complications  indicated 
by  body  wave  studies  (see  Fig.  22)  are  neglected, 
it  is  possible  to  explain  all  available  dis¬ 
persion  data  in  Central  Europe  by  slight  modifi¬ 
cations  of  model  STU-3, 

This  simplified  interpretation  was  justified 
until  it  became  clear  from  body-wave  results 
that  there  are  at  least  two  velocity  reversals 
in  the  upper  mantle,  probably  even  three.  In 
Fig,  28  (right)  the  phase  velocity  results  for 
the  line  Stuttgart-Oropa  are  plotted  as  a 
function  of  wave  period  (T)  or  wave  length  (L) , 
respectively.  The  latter  quantity  takes  it 
easier  to  visualize  the  corresponding  depth  of 
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penetration  for  these  signals.  Waves  with  the 
longest  periods  observed  will  "feel"  the  ve¬ 
locity  distribution  down  to  a  depth  of  about 
180  km.  The  fundamental  mode  and  particularly 
the  higher  mode  waves  should  therefore  be 
sensitive  to  velocity  inversions  within  the 
uppermost  part  of  the  mantle. 

Phase  velocities  for  two  models  which  are  based 
on  the  most  recent  travel-time  analysis  of 
explosion  and  earthquake  data  were  calculated 
and  are  shown  in  Fig,  28,  Model  DARMSTADT  003 
which  evolved  from  Model  ©  in  Fig.  15  pro¬ 
vides  a  satisfactory  explanation  of  the  dispersion 
measurements  for  the  period  range  from  20  to  120 
seconds.  It  should  be  kept  in  mind  that  the  phase 
velocities  for  lines  completely  outside  the  Alps 
are  on  an  average  higher  by  about  0.10  to  0.15 
km/sec,  thus  being  in  even  better  agreement  with 
the  theoretical  results. 

Model  SW-Europe  160  uses  the  velocity  distributions 
found  from  earthquake  travel  times  (see  Fig.  22) . 
Except  for  the  crustal  structure  which  in  the 
latter  case  has  beon  adopted  from  investigations 
on  the  Iberian  Peninsula  (PAYO,  1965)  there  are 
pronounced  differences  in  the  number  and  depth 
intervals  of  the  P  and  S  velocity  reversals.  It 
is  seen  in  Fig,  28  (right)  that  the  calculated 
phase  velocities  are  consistently  too  low  in  the 
intermediate  period  range.  This  is  partly  due  to 
the  different  crustal  structure  assumed*  The 
theoretical  results  should,  therefore,  not  be 
compared  to  data  from  Central  Europe.  In  addition 
the  mean  velocity  of  the  uppermost  100  km  seems 
to  be  too  low.  For  all  these  reasons  we  presently 
prefer  model  DARMSTADT  003, 
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c  iKm/sec) 


Fig,  28:  Comparison  of  experimental  Rayleigh  wave  disnoi"! 

results  with  theoretical  calculations  for  the  now 
models  DARMSTADT  003  and  SV7-Europe  160 
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4.3  Amplitudes  Attenuation  of  Surface  Waves 

(U.  WALTER) 

Since  amplitudes  of  body  waves  are  particularly 
difficult  to  interpret  because  of  source  factors, 
mode  conversion  at  interfaces,  scattering, 
spreading  losses  and  instrumental  effects,  a 
new  technique  has  recently  been  developed  for 
interpreting  the  frequency-dependent  amplitude 
decay  data  of  surface  waves  in  term3  of 
anelasticity  versus  depth  (ANDERSON  and 
ARCHAMDEAU ,  1964 >  ANDERSON,  BEN-MENAHEM  and 
ARCHAMBEAU,  1965). 

We  Fourier-analyzed  vertical  pendulum  records 
of  several  nuclear  explosions  in  the  atmosphere 
near  Novaya  Zemlya  and  also  of  a  number  of 
natural  earthquakes  in  order  to  measure  the 
phase  and  group  velocities  as  well  as  the 
rate  of  decay  of  Rayleigh  wave  energy  in  the 
spectral  band  of  about  10  to  100  seconds.  These 
data  are  used  to  determine  the  attenuation  in 
the  upper  mantle.  In  particular  the  eight  events 
listed  in  Table  4  were  studied i 


(1) 

23 

October 

1961 

(2) 

30 

October 

1961 

(3) 

5 

August 

1962 

(4) 

19 

September 

1962 

(5) 

25 

September 

1962 

(6) 

27 

September 

1962 

(7) 

22 

October 

1962 

(8) 

24 

December 

1962 

Sample  seismograms  for  the  second  event  as 
recorded  at  Copenhagen  (COP) ,  Stuttgart  (STU) 
and  Strasbourg  (STR)  are  reproduced  in  Fig.  29. 
Two  interfering  wave  trains  can  be  clearly 
discerned. 
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Fig,  29»  Sample  seismograms  for  tho  atmospheric  explosion 

near  Novaya  Zemlya  on  30  October  1961  as  recorded 
at  COP,  STU  and  STR 
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Out  of  those  eight  explosions,  the  third  and 
fourth  were  selected  for  a  more  extensive 
study.  The  long-period  vertical  records  of 
the  WWSS-Network  for  the  available  European 
stations  were  digitized  over  an  interval 
covering  the  entire  surface  wave  train.  Reading 
inaccuracies  could  not  bo  avoided  in  those 
parts  of  the  seismogram  where  the  amplitudes 
were  largest.  The  digitized  data  were 
harmonically  analyzed  and  the  calculated 
Fourier  amplitudes  wore  then  corrected  for 
different  peak  magnification  of  the  seismographs. 

In  Figs.  30  and  31  the  computed  amplitude  spectra 
for  selected  stations  are  depicted.  For  reference 
the  amplitude  response  of  the  30-100  system  then 
in  operation  at  the  V7WSS  stations  is  also  shown. 
Fairly  smooth  spectra  have  been  obtained  for 
epicentral  distances  up  to  abouu  2  500  km,  i.c. 
for  stations  located  on  the  Fennoscandian  shield. 
Severe  deformations  of  the  spectrum  are  observed 
as  soon  as  the  transition  occurs  between 
Fennoscandia  and  Central  Europe,  Changos  in 
crustal  structure  must-be  responsible  for  the 
observed  attenuation  at  the  short-period  end 
of  the  spectrum,  while  absorption  in  the  upper¬ 
most  mantle  is  the  cause  of  the  spectral 
degradation  for  periods  greater  than  about 
30  seconds. 

Following  TRYGGVASON  (1965)  an  attempt  has  been 
made  to  determine  the  dimensionless  quality 
factor  Q  for  Rayleigh  wave  propagation  across 
Europe.  In  the  period  range  of  around  40  to 
50  seconds  it  has  a  value  close  to  250,  A 
preliminary  interpretation  shows  that  Q  must 
increase  rapidly  as  the  period  decreases.  At 
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shorter  periods  and  wave  lengths  corresponding 
to  crustal  dimensions  the  attenuation  again 
increases/  but  the  scatter  of  the  data  is  too 
large  and  does  therefore  not  permit  any 
quantitative  conclusions* 

The  amplitude  decay  data  for  Love  (Q^)  and 
Rayleigh  (QR)  surface  waves  consistent  with 
available  body  wave  data  are  summarized  in 
Fig.  32*  Our  QR  average  for  Europe  (o)  is 
plotted  together  with  TRYGGVASON's  value  (x) 
on  top  of  the  data  points  by  ANDERSON, 
BEN-MENAHEM  and  ARCIIAMBEAU  (1965)  for  the 
paths  Toledo-Trinidad  and  Chile-Pasadena.  The 
specific  quality  factor  Q  is  shown  as  a 
function  of  period  (T) •  As  before  in  Fig*  28 
this  type  of  display  should  make  it  easier 
to  visualize  the  depth  of  penetration  if 
the  appropriate  fractions  of  wave  length 
are  considered.  Without  too  much  imagination 
there  is  clear  evidence  for  two  zones  of 
relatively  low  Q  in  the  upper  part  of  the 
mantle  both  for  Love  and  Rayleigh  waves, 
indicating  that  the  increased  attenuation  is 
primarily  due  to  losses  in  shear  (Qg)  . 

A  rough  estimate  shows  that  zones  of  low  Qs 
must  exist  at  depths  of  about  80  -  110  km 
and  150  -  220  km.  These  depth  intervals 
correspond  to  the  depths  found  for  the  two 
main  low-velocity  channels  discussed  above. 

We  conclude  from  these  results  that  the  simple 
Q  model  MM  8  proposed  by  ANDERSON,  BEN-MENAHEM 
and  ARCHAMBEAU  (1965)  is  not  adequate  to 
explain  the  observations  for  periods  smaller  than 
150  seconds.  The  Q  values  found  for  the  crust 
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32 »  Summary  of  attenuation  data  for  Love  (Q )  and 
Rayleigh  (Qn)  waves.  The  solid  dots  rcnrescnt 
data  by  ANDERSON,  BEN-MENAIIEM  and  ARCHAMBEAU  (1965) 
the  x  is  a  value  by  TRYGGVASON  (1965),  and  o  is 
our  average  value  for  Europe 
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are  much  lower  than  the  ones  deduced  from  the 
earlier  results  of  PRESS  (1964) »  Additional 
complice. c-ons  must  be  expected  for  the  depth 
range  immediately  below  the  crust-mantle 
transition  zone*  Much  more  shorter-period  data 
of  improved  quality  are  now  needed  to  answer 
all  these  unsolved  questions. 


5.  ANELASTIC  PROPERTIES  OF  THE  DEEPER  MANTLE 


5*1  Average  Q  of  Congressional  Waves  (Qp) 

(K.  FUCHS) 

A  rough  estimate  of  Qp  for  the  upper  mantle 
has  been  obtained  based  on  P  observations  of 
two  Hindu  Kush  intermediate  focus  earthquakes 
(focal  depth  h^OO  km)  as  recorded  at  European 
and  African  WWSS  network  stations.  The  ampli¬ 
tude  spectra  of  the  P  phases  for  the  two  shocks 
arranged  according  to  increasing  distance  are 
displayed  in  Fig.  33.  They  have  been  corrected 
for  instrumental  distortion.  Spectral  cor¬ 
rections  for  the  receiver  crust  have  not  yet 
been  applied  since  not  enough  reliable  in¬ 
formation  on  crustal  structure  was  available 
for  all  the  stations.  It  is  felt  that  the 
average  which  is  determined  from  a 
number  of  stations  with  a  wide  variety  of 
crustal  structures  should  not  be  too  much 
in  error  as  the  influence  of  the  crust  will 
cancel  out  to  a  first  approximation. 


If  it  is  assumed  that  the  intrinsic  Q  is 
independent  of  frequency  as  appears  to  be 
the  case  for  most  solid  materials  that  have 
been  measured  in  the  laboratory,  then  the 
amplitude  of  a  seismic  phase  may  be  expressed 
in  the  form* 

A  ( A,  h,to)  =  exp  (  -  T~  '  IF  ) 

^  Up 

In  this  equation  t  is  the  travel  time  along 
a  ray  starting  at  a  deep  focus  of  depth  h 
and  arriving  at  a  station  in  the  opi- 
central  distance  &  • 

Furthermore  if  it  is  assumed  that  the  fre¬ 
quency-dependent  amplitude  distortion  due  to 
anelasticity  is  about  the  same  for  all  rays 
above  the  depth  of  the  focus,  then  the  ob¬ 
served  differences  in  spectral  content  mast 
be  attributed  to  those  parts  of  the  ray  paths 
which  have  penetrated  into  the  deeper  mantle. 
More  precisely  only  the  propagation  along  rays 
through  an  earth  model  stripped  down  to  the 
level  of  the  focus  (h 200  km)  must  be  con¬ 
sidered.  To  this  end  the  stripped  earth  time 
t'  for  P  observed  at  a  station  with  opicentral 
distance  has  been  computed  for  different 
earth  models  using  the  same  ray  parameter  p 
in  the  stripped  earth  which  would  correspond 
to  the  opicentral  distance  A  and  the  focal 
depth  h  in  the  complete  model. 

With  the  above  assumptions  andco  -  Z7CV  the 
amplitude  spectrum  can  now  be  written  ast 

A(v)  -  exp  f-  Jj-tLy,  +  D(V) 

Qt 
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The  term  D  (  V  )  takes  care  of  the  source 
spectrum  which  presumably  is  flat  for  all 
rays  leaving  the  "focal  window"  within  a 
narrow  anglo  of  aperture  (8°  in  the  present 
case) .  D  (  V  )  also  includes  the  effect  of 
the  ray  segment  above  the  stripped  earth  at 
the  receiving  end. 


In  Fig,  33  it  i3  S2cn  that  the  logarithms 
of  the  normalized  amplitudes  show  a  roughly 
linear  dependence  on  frequency  V  within  the 
observed  spectral  range.  Therefore  to  a  first 
approximation  D  ( V  )  must  also  be  linearly 
dependent  on  frequency  within  that  range: 


D  (v)  s  -  E •  » 

The  measured  spectral  gradients  B  «  th  A  /v 
and  the  corresponding  stripped  earth  times  t' 
of  a  Joffreys-Bullen  earth  model  were  then 
used  to  determine  a  r.m.s.  best-fit  line  to 
the  equation: 


For  the  two  events  the  results  arc  a3  follows: 

6  July  1962:  =»  149  -  40  ,  E  a  1,6  -  2.4 

28  Jan,  1964:  ■»  201  -  86  ,  E  ■  1.5  -  3.0 

With  the  hic.i  error  limits  the  assumption  E  ^  0 
seems  to  be  justified,  i.e.  the  source  spectrum 
must  be  essentially  flat  within  the  frequency 
range  considered.  Again  fitting  a  r.m.s.  lino 
to  the  data  results  in: 

6  July  1962:  -  179  -  5 

28  Jan,  1964:  ■  256  -  11 
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It  follows  that  the  average  Q_,  value 

t  « 

(Qp®*220)  for  the  deeper  part  of  the  mantle, 
i.c.  for  depths  between  200  and  1600  km  Is 
much  lower  than  the  values  which  have  been 
reported  In  the  literature  so  far* 


5*2  Implications  on  Loss  Mechanisms 

A  complete  description  of  possible  loss 
mechanisms  In  the  earth  requires  the  measurement 
of  the  attenuation  of  compresslonal  and  shear 
waves.  According  to  KOVACH  and  ANDERSON  (1964) 
an  average  value  of  200  for  in  the  upper 
600  km  of  the  mantle  satifies  the  observational 
data.  The  average  of  600  in  the  whole  mantle 

and  the  value  for  of  200  in  the  upper  mantle 

leads  to  an  estimate  of  2200  for  the  average 
q]T  of  the  lower  mantle.  This  number  is  an 
order  of  magnitude  greater  than  the  0g  in  the 
upper  mantle.  The  observed  distribution  of  Qg 
with  depth  suggested  by  these  body  wave  results 
seemed  to  agree  with  the  Qg  distribution  inferred 
from  torsional  oscillation  data  and  Love  wave 
observations  of  ANDERSON  and  ARCIIAMBEAU  (1964). 

A  promising  approach  in  determining  the 
attenuation  of  compresslonal  waves  in  the  mantle 
has  been  suggested  by  KOVACH  (1967)  •  The 
logarithm  of  the  spectral  amplitude  ratio  of 
SKP  and  SKS  permits  the  attenuation  of  com- 
pressional  waves  (Qp)  in  the  mantle  to  be  com¬ 
pared  to  that  of  shear  waves  (Qg)  since  the 
effects  of  the  source  and  of  the  common  pro¬ 
pagation  paths  can  be  eliminated.  Preliminary 
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results  indicate  that  the  average  Qp  for  the 
whole  mantle  is  about  500,  i.o.  comparable 
to  the  average  of  600  in  the  mantle  and 
much  lower  than  has  been  assumed  until 
recently. 

Several  years  ago  CARPENTER  and  FLINN  (1965) 
presented  results  which  seemed  to  indicate 
that  the  following  proportionality  holds: 

Qp  °  ^ 

It  was  claimed  that  this  relation  with 
1.00<k<1.25  holds  for  P  signals  recorded 
at  teleseismic  distances  (A  > 20°)  from 
shallow  events.  Attenuating  effects 
of  both  the  source  and  the  receiver  system 
are  thus  includod  in  the  analysis.  The  data 
presented  in  this  report  suggest  that  except 
for  the  crustal  segment  of  the  ray  paths  the 
average  is  only  about  one  half  of  the  mean 
value  obtained  from  the  above  proportionality 
relation.  This  is  at  least  substantiated  by 
our  P  observations  in  the  period  range  2<T<20 
sec  and  for  epicentral  distances  out  to  60°, 
i.e,  down  to  depths  of  about  1500-1600  km. 

A  recent  study  by  TENG  (1966)  has  again  led 
to  the  generally  accepted  result  that  the 
lower  mantle  has  a  much  higher  Q  than  the 
upper  mantle.  But  the  rapid  increase  of  0p 
postulated  now  must  occur  at  a  depth  of 
about  950  to  1000  km,  considerably  deeper  than 
the  estimates  of  450  km  (ANDERSON,  BEN-MENAHEM 
and  ARCHAMBEAU,  1965)  and  650  km  (KNOPOFF, 

1964) »  By  assuming  no  losses  in  pure  compression 


Op  should  be  approximately  2,5  times  Qs,  thus 
ruling  out  major  nonadiabatic  or  thermal 
conduction  losses*  Under  this  assumption 
ANDERSON,  BEN-MENA1IEM  and  ARCHAMBEAU  (1565) 
have  derived  a  Qp  model  for  the  mantle  from 
their  Qs  structure  that  is  obtained  by 
inverting  surface-wave  data. 

In  agreement  with  TENG  (1968)  wo  find  that  our 
results  are  closer  to  Qp«^Qs,  and  that  con¬ 
sequently  losses  in  pure  compression  are  not 
negligible.  Similar  conclusions  have  also  been 
reached  by  KANAMORI  (1967  a+b)  for  periods  of 
1-2  seconds.  Even  if  allowance  is  made  for 
experimental  uncertainties,  the  discrepancy 
between  the  body-wave  and  surface-wave  results 
seems  to  be  significant.  It  is  thus  indicated 
that  around  1-2  sec  period  may  be  smaller 
by  a  factor  of  2  compared  to  around  25-50  sec 
period,  i.e.  that  decreases  with  increasing 
frequency. 

These  results  would  suggest  that  the  approximation 
of  a  frequency-independent  Q  does  not  hold  over 
a  wider  frequency  band.  A  frequency  dependence 
of  Q  cannot  be  established  at  this  point  because 
of  the  paucity  of  data.  But  if  the  result 
QpA'Qg  is  accepted,  losses  in  pure  compression 
are  no  longer  small  or  negligible  compared  to 
losses  in  shear,  and  loss  mechanisms  such  as  for 
example  thermoelasticity  cannot  be  ruled  out 
completely. 
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7.  CHRONOLOGICAL  BIBLIOGRAPHY 

The  following  reports  on  work  performed  during  the  dura¬ 
tion  of  this  contract  have  been  completed  and  delivered 
to  the  European  Office  of  Aerospace  Research  (EOAR) x 


Tyne  of  Report 

Period 

Covered 

Letter  Planning  Report 

(Administrative  Report) 

No.  1 

1 

May  1965 

- 

31 

May 

1965 

Quarterly  Narrative  Progress  Report 
(Scientific  Report)  No,  1  1 

May  1965 

- 

31 

July 

1965 

(1  H 

No,  2 

1 

Aug.  1965 

- 

31 

Oct. 

1965 

II  If 

No.  3 

1 

Nov.  1965 

- 

31 

Jan, 

1966 

Annual  Summary  Report 

No.  1 

1 

May  1965 

- 

30 

April 

1966 

Quarterly  Narrative  Progress  Report 
(Scientific  Report)  No,  4  1 

May  1966 

- 

31 

July 

1966 

If  II 

No.  5 

1 

Aug.  1966 

- 

31 

Oct. 

1966 

•1  If 

NO.  6 

1 

Nov,  1966 

- 

31 

Jan  • 

1967 

Progress  Report 

NO.  7 

1 

Febr.1967 

30 

April 

1967 

91  IV 

No.  8 

1 

May  1967 

- 

31 

July 

1967 

II  II 

NO.  9 

1 

Aug,  1967 

- 

31 

Oct. 

1967 

Final  Scientific  Report 

1 

May  1965 

- 

30 

April 

1968 

These  reports  were  supplemented  by 
Special  Scientific  Report  No,  1,  dated  20  May  1966  ,  on 

"  Seismic  Ground  Noise  and  Wind  at  the 
Seismological  Observatory  Graefenborg  " 
by  H.  BERCKHEMER  and  B.  AKASCHE  . 

As  a  result  of  this  contract  work  a  number  of  publications 
will  soon  be  available*  The  first  paper  has  }ust  been 
published! 

DOHR,  G,  and  K,  FUCHS,  Statistical  evaluation  of  deep 
crustal  reflections  in  Western  Germany, 
Geophysics,  32,  951-967,  1967 
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Throe  papers  are  about  to  be  published: 

FUCHS,  K,,  Das  Reflexions-  und  Transmissionsvcrmogen 
eines  gcschichteten  Mediums  bei  bclicbi- 
ger  Tiefenvorteilung  dor  elastischen  Mo- 
duln  und  der  Dichte  f  Jr  schriigen  Einfall 
ebener  Wellen, 

Ztschr. f .Goophysik,  34  (in  press) ,  1969 

MUELLER,  St.,  J.  ANSORGE,  D.  MAYER-ROSA  and  K.  FUCHS, 

Evidence  for  several  velocity  reversals 
in  the  upper  mantle, 

Ztschr. f .Goophysik,  34  (in  press),  1968 

FUCHS,  K,,  On  the  properties  of  deep  crustal  re¬ 
flectors, 

Ztschr. f . Goophysik,  34  (Jin  press),  1969 

Four  more  papers  with  the  following  titles  are  in  pre¬ 
paration: 

(1)  "  Structure  of  the  crust  and  uppermost  mantle 

in  the  eastern  part  of  France  " 

(2)  "  Das  DurchlaQspektrum  der  Erdkruste  in  Europa 

fUr  langperiodi3che  seismische  Raumwellen  " 

(3)  "  Energieverluste  3eismischer  Raumwellen  in  der 

Erdkruste  und  im  obersten  Erdmantel  " 

(4)  "  Attenuation  of  compressional  waves  in  the 

earth's  mantle  " 
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(%.  MIT  ■  ACT  " 

The  attenuation  of  seismic  v;aves  in  the  crust  and  upper  mantle  has  bee 
studied  by  various  techniques  in  order  to  determine  the  specific 
quality  factor  Q,  Andlastic  deficiencies  qf  the  propagation  medium 
can  only  be  separated  after  corrections  are  made  for  the  effects  of 
instrumental  distortion,  of  the  "receiver . crust"  and  the  "source  crustf 
the  lossless  propagation  through  the  mantle  and  the  source  mechanism,  ( 

Seismic  refraction  data  have  boon  compiled  for  a  number  of  "source"  and 
"receiver"  regions  in  Europe,  A  refined  ve lo city -dentil  distribution 
with  two  velocity  reversals  and  an  average  value  of  Q  25 670  have  boon 
obtained  for  the  crust  of  the  Bohemian  Massif  in  Central  Europe,  The 
Q_j  of  the  uppermost  mantle  was  found  to  be  about  250,  Crustal  transfer 
functions  were  calculated  to  tost  models  of  crustal  structure  derived 
by  other  methods.  The  radiation  patterns  for  two  typos  of  explosive 
sources  have  boon  determined.  Their  shape  is  strongly  affected  by  the 
source  mechanism  and  crustal  structure, 

A  new  velocity-depth  structure  of  the  upper  mantle  in  Europe  has  been 
derived  from  travel-time  observations  of  large  explosions  and  earth¬ 
quakes,  At  least  three  velocity  reversals  have  been  found  down  to 
depths  of  about  220  km.  They  correlate  well  with  the  distribution  of 
earthquake  focal  depths,  Surface-wave  dispersion  measurements  are  also 
in  agreement  with  those  now  findings.  Amplitude  decay  data  for  surface 
waves  indicate  that  .zones  of  low  Qr  must  exist  at  depths  whe re  the 
low-vclocity  channels  wore  found,  0 
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13 ,  Abstract  (continued) 


Finally  an  average  value  of  Qp«s  220  for 
the  deeper  mantle  has  been  obtained  based 
on  observations  of  two  intermediate  focus 
earthquakes,  A  comparison  with  corresponding 
data  shows  that  QpwQs,  i.o.  that  losses 
in  puro  compression  cannot  be  neglected. 
Alternatively  it  is  suggested  that  Q  may 
bo  frcquer.cy-dopendent,  at  least  over  an 
extended  frequency  band. 
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